Gate-tunable bandgap in bilayer graphene 
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The tight-binding model of bilayer graphene is used to find the gap between the conduction and 
valence bands, as a function of both the gate voltage and as the doping by donors or acceptors. The 
total Hartree energy is minimized and the equation for the gap is obtained. This equation for the 
ratio of the gap to the chemical potential is determined only by the screening constant. Thus the 
gap is strictly proportional to the gate voltage or the carrier concentration in the absence of donors 
or acceptors. In the opposite case, where the donors or acceptors are present, the gap demonstrates 
the asymmetrical behavior on the electron and hole sides of the gate bias. A comparison with 
experimental data obtained by Kuzmenko et al demonstrates the good agreement. 

PACS numbers: 73.20.At, 73.21.Ac, 73.43.-f, 81.05.Uw 



INTRODUCTION 

Bilayer graphene has attracted much interest partly 
due to the opening of a tunable gap in its electronic 
spectrum by an external electrostatic field. Such a phe- 
nomenon was predicted in Refs. [H, Q and can be ob- 
served in optical studies controlled by appl yin g a gate 
bias d, II i, i, S i, S, In Refs. P, [l^, within 

the self-consistent Hartree approximation, the gap was 
derived as a near-linear function of the carrier concen- 
tration injected in the bilayer by the gate bias. Recently, 
this problem was numerically considered in Ref. [13| us- 



ing the density functional theory (DFT) and including 
the external charge doping involved with impurities. The 
DFT calculation gives the gap which is roughly half of 
the gap obtained in the Hartree approximation. This dis- 
agreement was explained in Ref. [l3| as a result of both 
the inter- and intralayer correlations. 

In this work, we study this problem within the same 
Hartree approximation as in Refs. 11 , l2| , but including 
the effect of external doping. We calculate the carrier 
concentration on both sides of the bilayer considering the 
case, where the carrier concentration in the bilayer is less 
than 10^^ cm~^. Then, we minimize the total energy of 
the system and find self-consistently both the chemical 
potential and the gap induced by the gate bias. Our 



11 



results completely differ from those obtained in Refs 
12|, where the external doping is disregarded. In the 
presence of dopants, the dependence of the gap on the 
carrier concentration, i.e., on the gate voltage, exhibits 
an asymmetry at the electron and hole sides of the gate 
bias. 



TIGHT-BINDING MODEL OF BILAYER 
GRAPHENE 



by solid lines, and in the other layer, e.g., ai and 6i, 
by the dashed lines. The atom a (ai) differs from b (bi) 
because it has a neighbor just below in the adjacent layer, 
whereas the atom b (bi) does not. 

Let us recall the main results of the Slonchewski- 



Weiss-McClure model [ij, |15|. In the tight-binding 
model, the Bloch functions of the bilayer are written as 



j 

= ^ I] ^'""'Mi^j + a - r) 



(1) 
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where the sums are taken over the lattice vectors a.j and 
N is the number of unit cells. Vectors a and c connect 
the nearest atoms in the layer and in the neighbor layers, 
correspondingly. 




FIG. 1: Bilayer lattice 



The graphene bilayer lattice is shown in Fig. [TJ Atoms 
in one layer, i. e., a and b in the unit cell, are connected 



For the nearest neighbors, the effective Hamiltonian in 



the space of the functions ([T]) can be written as 



H{k) 



fU + A 70/* 71 74/ \ 

70/ U-A 74/ 73/* 

71 74/* -U + A 70/ 

V 74/* 73/ 70/* -U-AJ 



(2) 



where / = 70 [e*'^'-'' + 2e-*'=-''/2 cos (fcj^a\/3/2)] . The val- 
ues of hopping integrals 70,71,73,74, and A are given in 
the Table I. The largest of them, 70, determines the band 

TABLE I; Parameter values of the electron spectrum in eV. 



Parameter 


Experiment [10] 


DFT calculation [16] 


70 


3.16 ±0.3 


2.598 ±0.015 


71 


0.381 ± 0.003 


0.34 ± 0.02 


73 


0.38 ± 0.06 


0.32 ± 0.02 


74 


0.14 ±0.03 


0.177 ±0.025 


A 


0.022 ± 0.003 


0.024 ±0.01 



dispersion near the K point in the Brillouin zone where 
the matrix element 70/ can be expanded as 

70/ = v{ikx - ky), 

with a velocity parameter v — 8700/2. The parameters 
73 and 74 giving a correction to the dispersion are less 
than 7o by a factor of 10. The parameters 71 and A re- 
sult in the position of levels at K, but A is much less 
than 71. There is in addition the parameter U induced 
by the gate voltage and associated with the asymmetry 
of two layers in the external electrostatic field. This pa- 
rameter presents the potential energy —edE between two 
layers, where d is the interlayer distance and E is electric 
field induced both by the gate voltage and the external 
dopants in the bilayer. 

The parameter U as well as the chemical potential /i 
should be self-consistently calculated for the given gate 
voltage. For this purpose, we can held only the parame- 
ters 7o and 71, neglecting the small effect of 73, 74, and 
A on the gap U . In this approximation, the effective 
Hamiltonian can be written in the simple form 



i/(k) 



f U vk+ -fi \ 
vk- U 

71 —U ffc_ 



V 



vk4 



(3) 



Uj 



where k± = ^ik^ — ky in the vicinity of the K points. 
The Hamiltonian gives four energy bands: 



e^Aq) = ± ( f - 



W 



-W 



1/2 



1/2 



(4) 




FIG. 2: Band structure of bilayer. 



where 



1/2 



H The min- 



and we denote q"^ = {vkY . 

The band structure is shown in Fig. 
imal value of the upper energy Si is yJU"^ + 7i • The 
£2 band takes the maximal value \U\ a.t q — Q and the 
minimal value U = 7i|C/|/\/7i + 4C/2 at q^ "^^V^ill + 
2t^^)/(7i + 4C/^). Because the value of U is much less 
than 71, the distinction between U and U is small and 
the gap between the bands £2 and £3 takes approximately 
the value 2\U\. 



EIGENFUNCTIONS AND CARRIER 
CONCENTRATION 

The four eigenfunctions C corresponding with eigen- 
values (UJ of the Hamiltonian ([3]) read 



/ ((7 -£„)[(£„ + C/)2-g2]N^ 

2 „2l 



c = 



1 
c 



-q-[{en + UY-q 



\ 



7i9+(C/-£„) 
where the C norm squared is 



(5) 



[(e„ + Ur~q'r[{en-Uy + q'] 
+jl{sn-Uf[{en + Uf +q^]. 

As seen from Eqs. ([T]), the probability pi to find an 
electron, for instance, on the first layer is 



Pi = \Ci 



\' + \c,\\ 



where the subscript Ci numerates the elements of the 
column dS]). 

We assume, that carriers occupy only the bands £2,3, 
so the chemical potential n and the gap 2\U\ are less 
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than the distance between the bands £i and £2, i-e., 
(|^|,2|C/|) < 71. The electron dispersion for the £2,3 
bands can be expanded in powers of q^: 



£^(g) = C/2-4^g2 + %, 
7f It 



(6) 



where n = 2 stands for the electron conductivity and 
7T, = 3 for the hole conductivity. Then, for ^ 4C/^, 
we can omit here the second term and use the simple 
relations 

q^=^,[el-U'Yl' (7) 

neglecting the small effect of the " mexican hat" . 

Keeping only the leading terms, one can find with the 
help of Eq. (O, that the probabilities pi,2 to find an 
electron on the layers are proportional 

Pi = |Cip + |C2pocg*^ = 7f(£2-[/2)3/2, 

P2 = |C3p + |C4p«g^?(C/-£„)' 

Therefore, the normalized probability to find an electron, 
for instance, on the first layer can be written as 



Pi 



(£2-C/2)3/2 



(£2 -C/2)3/2 + (e2 _[/2)l/2(t/_£„) 

= (£„ + C/)/2£„ . (8) 



Within the approximation ([7])- ([5]), many observable ef- 
fects can be analytically evaluated for the intermediate 
carrier concentration, AU"^ <^ 71 v^^^ _ {72 ^ ^2_ 

At zero temperature, the carrier concentration on the 
sides of the bilayer is found with the help of Eq. ([5]) as 



ni.2 



Pi, 2 qdq 



npU 
271 



[\/x^~l ± In {x + \/x^ - 1] 



(9) 



where the limits of integration arc q = and the chemical 
potential /x and we set 



no = -ff/TTh^v^ = 1.03 X 10"cm~2 , 

x = |m/C/|. 



(10) 



For the total carrier concentration n in the bilayer, we 
obtain 



fc2 2 

nn V 



71 



1 



(11) 



MINIMIZATION OF THE TOTAL ENERGY 

In order to find the chemical potential /i and the gap 
2\U\ at the given gate voltage 



-edE — edwEr, 



(12) 



^1 



E 

- 

d 



-n 



w 
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FIG. 3: Electrostatic model; d is the interlayer distance, 
is the wafer thickness. 



we minimize the total energy containing both the energy 
y^'^) of the carriers and the energy V^^"^ of the electro- 
static field. Within the Hartree approximation, when 
no electron correlations are taken into account, the filled 
bands do not contribute into the energy of the system, 
because the electron charge of the filled bands is com- 
pensated by the ion charge and this energy have to be 
considered as the ground state energy. The excitation 
energy is owes its origin to the carriers in unfilled bands. 
Electrons in the £2 band or holes in the £3 band con- 
tribute in the total energy of the system the energy 



\eniq)\qdq 



271 



[x\/x^~l + In (a: + V x'^ - 1)] . 



The energy of the electrostatic field (see Fig. [3]) 

= ±^dE^ + e^d^El) 
an 



(13) 



(14) 



can be written in terms of the carrier concentrations with 
the help of relations 

4ne{ni-Ni)=E and 4Tre{n - N) ^ SyjEyj , (15) 

where is the dielectric constant of the wafer, Ni and 
A^2 are concentrations of the acceptor or donor impurities 
on the left and right layers, correspondingly, whereas the 
total dopant concentration on the bilayer is N = Ni + 
All these numbers are supposed to be positive or 
negative for the electron or hole doping correspondingly. 
Let us emphasize that the dielectric constant e of bilayer 
graphene depends on the substrate. For simplicity, we 
put e = 1 in the definition (fTl]) . 

We seek the minimum of the total energy 

+ + xi^eVg + edE + ed^uE^) 

as a function of U, /x, and the Lagrange multiplier A. 
Differentiation with respect to A gives the gate voltage 
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constraint ((12)) . Minimization with respect to U and x 
gives 

47re^[ni - Ni)niud+ {n~ N)nudw/^w 
+yi'=) + Ai:e^\{niud + n„d„/e„) = 

and the similar equation with a substitution u x, 
where the subscripts u and x note the derivatives with 
respect to the corresponding variables. The Lagrange 
multiplier A can be excluded from these two equations. 
Then, the equation obtained should be expanded in d/d-^j, 
since the thickness d of the bilayer is much less than the 
thickness dw of the substrate. 

Thus, we obtain the following equation: 



V "a; nu J 



(16) 



Let us emphasize, that this equation is invariant under 
the simultaneous sign change in ni^2 and Ni^2, that ex- 
presses the charge invariance of the problem. At the fixed 
external doping iVi^2, the gap on the electron and hole 
sides of the gate bias is not symmetrical. 

The derivatives in Eq. are calculated with the 

help of Eqs. ©-((Tni). As a result, Eq. HI]) becomes 



, 71^2 

' Una 



y/x^-l± |/(a:) + 



xf{x) 



A[xf{x) - VxP 



17) 



with the function /(x) = In (a; + \/x^ — 1) and the di- 
mensionless screening constant 



A: 



e'^jid 
Jhvf 



(18) 



For the parameters of graphene d = 3.35 A , 71 = 0.381 
eV, and v = 1.02 x 10* cm/s, we get A = 0.41. 



THE GAP IN UNDOPED AND DOPED BASED 
BILAYER 



(i) undoped bilayer 
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FIG. 4: The gap in units of 71 = 0.381 eV versus the carrier 
concentration in absence of doping (dashed line) and for the 
hole doping level = —2.5 x 10^^ cvcT^ (solid line); the pos- 
itive (negative) values of n correspond to the electron (hole) 
conductivity. The difference between values of n marked as 
"gate bias = 0" and "minimal dc conductivity" is 2A'2. 



dependence is shown in Fig. 3] in dashed lines, it is sym- 
metrical on the electron and hole sides. 



With the help of Eq. ([TOl) . we obtain the chemical 
potential as a linear function of the carrier concentration 



7ia^o 

/i = , ^ n . 



- 1 



(20) 



where n is positive (negative) for the electron (hole) con- 
ductivity. 

We can compare Eq. (fTS)) with the corresponding re- 
sult of Ref. iUl]: 



2|f//n| 



2^ 



1 + 2AM+Aln;^ 

no \n\ 



-\ -1 



(21) 



First, let us consider an ideal undoped bilayer with 
Ni — N2 — 0. We get a nonzero solution for [/, if the 
right-hand side of Eq. ([T7|) vanishes. This condition is 
fulfilled only for the sign "-" in Eq. (fTT)) . that defines the 
polarity of the layers [see Eq. We obtain solution 

X = xq = 6.61. According to Eq. pT|) . the gap as a 
function of the carrier concentration takes a very simple 
form: 



2\U/n\ 



271 



no\Jxl - 1 



1.13 X 10""meV • cm^ , (19) 



where the right-hand side does not depend at all on the 
gate bias, but only on the screening constant A. This 



Both equations give approximately the same results at 
\n\ ~ O.Itt-o ~ 10^^ cm^^. However, contrary to Eq. 
(HH]), Eq. (pij) contains the carrier concentration in the 
right-hand side giving rise to the more rapid increase in 
the gap with |7i| <C no- This increase also contradicts to 
the DFT calculations 

Two reasons can result in the disagreement of our the- 
ory with Ref. ll|. First, in Ref. ll|, the filled bands are 
supposed to create the electric field in the bilayer, that 
is incorrect as explained in the previous section. Second, 
the minimization should be done with respect two vari- 
ables and t/, however, only one of them seems to be 
exploited in Ref. 11 1. 
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(ii) doped bilayer 

For the bilayer with the acceptor or donor dopants, Eq. 
p?| presents a solution w = ^'^iNijUfiQ as a function of 
X. We obtain, evidently, the small values of w for x close 
to a;o = 6.61. Since xq ^ 1, we can expand the function 
in the right-hand of Eq. \YI\ in Xjx. In this region of the 
relatively large |[/|, we find again with the help of Eqs. 
PT|) and pT|) the linear dependence 

2ti 

2|J7| = |n-2iV2|^^ (22) 
= 1.13|n-27V2| X 10"" meV-cm^. 

The value of the carrier concentration n = 2N2 corre- 
sponds to the zero bias voltage, where U — (see Fig. 
E]). Therefore, in contrast with the undoped case, the gap 
demonstrates the asymmetrical behavior on the electron 
and hole sides. If the bilayer contains acceptors with con- 
centration N2, the gap decreases linearly with the hole 
concentration and vanishes, when the gate bias is not 
applied and the hole concentration equals 2N2. Starting 
from this point, the gap increases and, thereafter, be- 
comes again small (equals zero in Fig. 2]) at the carrier 
concentration corresponding to the minimal value of the 
dc conductivity, where n — 0. Therefore, the difference 
(1.56 X lO^^cm"^ in Fig. |4|) observed in Refs. and 
[10| between these two values of carrier concentrations, 
at the zero bias and at the minimal conductivity, gives 
directly the donor/acceptor concentration {2N2) on the 
layer close to the substrate. Then, for the gate bias ap- 
plied in order to increase the electron concentration, the 
gap is rapidly opening with the electron appearance. 

We see, that the asymmetry arises between the electron 
and hole sides of the gate bias. This asymmetry can 
simulate a result of the hopping integral A in the electron 
spectrum [isj . In order to obtain the gap dependence 
for the case of electron doping, iV2 > 0, the reflection 
transformation n — > —n has to be made. This case is 
shown in Fig. [5] where the experimental data from Ref. 
[lot are displayed. 

The gap in the vicinity of the minimal conductivity 
value could reach indeed a finite value due to several rea- 
sons. One of them is the form of the "mexican hat" shown 
in Fig. [2] Second, the trigonal warping is substantial at 
low carrier concentrations. Finally, the graphene electron 
spectrum is unstable with respect to the Coulomb inter- 
action at the low momentum values. For the graphene 
monolayer, as shown in Ref. [l^ . the logarithmic cor- 
rections appear at the small momentum. In the case of 
the bilayer, the electron self-energy contains the linear 
corrections, as can be found using the perturbation the- 
ory. The similar linear terms resulting in a nematic order 
were also obtained in the framework of the renormaliza- 
tion group [20j . 



minimal 




n ( 10^^cm"2) 

FIG. 5: The gap in eV versus the carrier concentration for the 
electron doping with the concentration N2 = 0.78 x 10^^ cm~^ 
(our theory); the positive (negative) values of n correspond 
to the electron (hole) conductivity; squares are experimental 
data [13]. 



CONCLUSION 

The gap opening in the gated graphene bilayer has an 
intriguing behavior as a function of carrier concentration. 
In the presence of the external doping charge, i.e. donors 
or acceptors, this function is asymmetric on the hole and 
electron sides of the gate bias and it is the linear function 
only for the large gate bias. A difference between two 
values of carrier concentrations, i. e. at the zero bias and 
at the minimal conductivity, gives directly the sign and 
concentration of the charged dopants on the bilayer. 
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